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Abstract 

Data from J/ -0 and ip' production m. p — A collisions are used 
to determine the cross section for absorption of pre-resonance char- 
monium in nuclear matter. The J/ij: suppression in O — Cu, O — U 
and S — U collisions is fully reproduced by the corresponding nuclear 
absorption, while Pb — Pb collisions show an additional suppression 
increasing with centrality. We study the onset of this change in terms 
of hadronic comover interactions and conclude that so far no con- 
ventional hadronic description can consistently account for all data. 
Deconfinement, starting at a critical point determined by central S—U 
collisions, is in accord with the observed suppression pattern. 
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1. Introduction 

The suppression of J /ifj production in nucleus- nucleus collisions was pro- 
posed ten years ago as a signal of colour deconfinement [T] . When a suppres- 
sion was observed m. O — U and later m. S — U collisions it was noted 
that already p — A interactions result in reduced J/ip production |Sj. Follow- 
ing recent studies of charmonium production in hadron-hadron interactions 
lUEj, J/'ip suppression 'm. p — A collisions can be understood as absorption 
of a pre- resonance cc — g state in nuclear matter The onset of colour de- 
confinement must then lead to an additional suppression of J/ip production, 
beyond the mentioned nuclear absorption; in other words, it must result in 
a reduction beyond what is already observed in p — A interactions [Zj. 

The "anomalous" J/ip suppression recently reported by the NA50 col- 
laboration [H-fni therefore triggered considerable excitement and already 
quite a number of tentative explanations [li|-[17j. While we fully share this 
excitement, we beheve that present data are sufficient to carry out a quan- 
titative analysis of charmonium suppression in nuclear collisions. There now 
exist quite extensive results on J/ ip production in p — A interactions for the 
determination of nuclear absorption, and the introduction of a zero degree 
calorimeter in the experimental set-up of NA50 specifies the underlying col- 
lision geometry much better than before. Combining the two in a systematic 
analysis leads to some clear-cut conclusions and also shows what further data 
is necessary to corroborate these. 

We begin our analysis, in Section 2, with a study of pre-resonance char- 
monium absorption in p — A collisions, using the Glauber theory of nuclear 
interactions. Apart from the well-established nuclear profile, this formalism 
contains as only parameter the absorption cross section aabs in nuclear mat- 
ter, which can thus be determined with a precision limited only by that of 
the p — A data. Turning to the collision geometry of nucleus-nucleus in- 
teractions, we establish in Section 3 the relation between the centrality of 
the collision, the associated transverse energy and the associated number of 
spectator nucleons. 

Using these results, we show in Section 4 that pre-resonance absorption 
fully accounts for the J/ip suppression observed in O — Cu, O — U and S — U 
collisions. In contrast, Pb — Pb collisions show considerable J/ip suppression 
beyond nuclear absorption. For the tp', there is increased absorption already 
in S — U collisions. 

In Section 5, we then investigate if the observed J/ip and ip' suppression 
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can be understood in "conventional" terms, i.e., as dissociation in a system 
of hadronic comovers produced in the course of the coUision. We find that 
the anomalous J / ip suppression in Pb — Pb and only nuclear absorption in 
S — U collisions make this impossible. In Section 6, we consider charmonium 
suppression by local colour deconfinement, assumed to begin at a critical 
density just above that of central S — U collisions. The resulting suppression 
pattern agrees with the Pb — Pb data. 

2. Nuclear Absorption in p — A Collisions 

Charmonium production in hadron-hadron collisions proceeds through 
parton fusion, at high energies predominantly gluon fusion, to form a cc pair 
|18j . Because of the high mass rric of the charm quark, this process occurs 
almost instantaneously, with a formation time Tc — (2mc)~^ — 0.07 fm in the 
cc rest frame. The cc pair is generally in a coloured state; to neutralize its 
colour and form a singlet cc state J/ijj or ip' requires a considerably longer 
time; virtuality estimates [Hj give about 0.3 fm, again in the rest frame of the 
cc. The analysis of recent experiments ^H] has shown the role of higher Fock 
space components in charmonium hadroproduction jH Ej. For production 
at low transverse momentum this suggests that the pre-resonance state of 
charmonia in the first 0.3 fm is a colour octet cc coupled with a soft collinear 
gluon to neutralise the overall colour of this cc—g system [H]. In the presently 
accessible kinematic region of J/ip production hj p — A collisions [xp > 0), 
the target nucleus sees only the passage of the pre-resonance state; physical 
charmonium states are formed outside the nucleus. 

The suppression of J/ip production in p — A collisions should thus be 
understood as pre-resonance absorption in normal nuclear matter. The size 
of the pre-resonance state is determined by the charm quark mass and the 
confinement scale; hence it is the same for J/ip and ip'- This accounts nat- 
urally for the equal suppression observed for the two states, which would 
be impossible for physical resonances whose transverse areas differ by more 
than a factor three. Theoretical estimates suggest for the absorption cross 
section of the cc — g on nucleons aabs — 6 — 7 mb We want to see to what 
extent it is possible to determine this cross section in a systematic study of 
presently available p — A data. 

In Glauber theory, the survival probability for a. J /ip produced ina.p — A 
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collision is given by 

= J (fb dz pA{b,z)exp l^-{A- 1) dz' pAib, z') aabs^ ■ (1) 

Here pA is the nuclear density distribution, for which we take the standard 
three-parameter Woods-Saxon form with parameters as tabulated [201; it is 
normalized to unity, with 

J d% dz pA{b, z) = 1. (2) 

In Eq. (^, the integration runs over the impact parameter b of the incident 
proton and the production point z of the pre-resonance state; the integral 
in the weight function covers the path of this state. The suppression is thus 
fully determined by the absorption cross section aabs in nuclear matter. 
Experimentally, the J/ip survival probability 

qexp _ exp I Aexp /o\ 
•^pA = '-'pA I ^'^pN ! 

is determined by comparing p — A with p — p or p — D results. We start 
with the latter, since here there are data at 200 [NA38], 450 [NA38/51] and 
800 GeV/c [E772]; for a recent compilation, see (TU]. In Table la, we hst the 
corresponding survival probabities. The errors of the 800 GeV/c data [21] are 
significantly smaller than those at the other energies because here p — A and 
p — D collisions were studied in the same experiment. Using the 450 GeV/c 
NA51 p — D cross section |22] to normalize the 450 GeV/c NA38 p — A data 
|2S1 results in larger errors. To normalize the 200 GeV/c p — A data, the 
450 GeV/c p — D cross section moreover has to be rescaled JU]; using the 
experimentally known energy dependence of J/ip production. - We now use 
the A-dependence of these results to determine the nuclear absorption cross 
section for pre-resonance charmonium. 

For the combined data in Table la we obtain the best fit with 

aabs = 7.3 ± 0.6 mb; (4) 

it gives a x^/d-f- — 1.4, with the error corresponding to a confidence level 
of 95% (see Fig. 1). Analyzing the 800, 450 and 200 GeV/c data separately 
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gives clear minima at 7.3, 7.4 and 7.1 mb, respectively. The results thus 
are fully compatible; they indicate that the A-dependence does not vary 
with incident energy in the presently studied range. This will presumably 
change when the beam energies become much higher (RHIC and LHC), since 
then nuclear shadowing of the gluon distributions is expected to enhance the 
suppression j^l- The survival probabilities obtained with the cross section 

are included in Table la. In Fig. 2, we compare the experimental and the 
nuclear absorption results; the agreement is very good. 

We have also carried out the corresponding analysis of the 450 and 200 
GeV/c p — A results with the NA51 p — p data as reference. The resulting 
o'abs = 6.4 ± 0.6 mb is somewhat lower, but still compatible with Eq. (j3J. 
However, here the large errors of the data make it difficult to define limits, 
so that the quoted range in this case corresponds to twice the minimum x^- 

So far, the nuclear suppression cross section was obtained with p — D 
or p — p data as reference. There is an alternative method which does not 
require these p—N values. The nuclear target ratios SpA/ SpB depend directly 
on aabs- We will here compare the results on different targets to those on 
carbon, since for this case there exist also 450 GeV/c data with relatively 
small errors The ratios at 800 and 450 GeV/c are listed in Table lb; 
they provide the best fit for aabs = 7.0 ±1.5, with x^/d-f- — 0.9 and the error 
determined by a 95% confidence level. The distribution is included in Fig. 
1; again the different experiments also have their individual minima at 7 
mb. We thus find that the value of aabs obtained through nuclear ratios is 
in good agreement with that obtained from Table la. Since the value Q is 
the most precise, we will use it in the remainder of our analysis. 

Next we turn to a comparison of J/ip and ip' production in p—A collisions. 
If the suppression observed in these reactions is indeed due to pre-resonance 
absorption, the ratio ip'/ {J/'^) should remain for all A at its value in nucleon- 
nucleon interactions, where it was found to be (1.81 ± 0.38) x 10~^ and 
independent of the c.m.s. collision energy from 20 to almost 2000 GeV fIE\ . 
Fitting the p — A results compiled in jTUj by the form 

B{J/'iIj ^ H+^i-) apA^j/^ 

we obtain (a^/ — a^) = 0.00 ±0.02 with a 95% confidence level; this rules out 
variations of more than 10% between p — p and p — U collisions. The value 
c = (1.74 ±0.07) X 10^^ agrees very well with the mentioned nucleon-nucleon 
average. 
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p 


nexp 


Sqi 




[GeV/c] 


(7.3 ±0.6 mb) 


p-C 


800 


0.851 ± 0.013 


0.867 ± 0.010 


p-Ca 


800 


0.806 ± 0.009 


0.784 ± 0.014 


p-Fe 


800 


0.756 ± 0.010 


0.753 ± 0.015 


p-W 


800 


0.619 ± 0.013 


0.650 ± 0.021 


p-C 


450 


0.85 ± 0.10 


0.87 ± 0.01 


p-Al 


450 


0.76 ± 0.09 


0.81 ± 0.02 


p-Cu 


450 


0.75 ± 0.08 


0.74 ± 0.02 


p-W 


450 


0.67 ± 0.08 


0.65 ± 0.02 


p-Cu 


200 


0.75 ± 0.17 


0.74 ± 0.02 


p-W 


200 


0.65 ± 0.05 


0.65 ± 0.02 


p-U 


200 


0.63 ± 0.13 


0.61 ± 0.02 



(a) J/ip suppression bom pA/pD [2,10,22] 





Pbeam 

[GeV/c] 


Qcxp 1 nexp 
'-'pA / '-'pC 


npA 1 npC 
'-'Gl 1 '-'Gl 

(7.3 ±0.6 mb) 


Ca/C 


800 


0.947 ± 0.024 


0.904 ± 0.007 


Fe/C 


800 


0.888 ± 0.025 


0.869 ± 0.009 


W/C 


800 


0.727 ± 0.026 


0.750 ± 0.015 


Al/C 


450 


0.925 ± 0.055 


0.931 ± 0.005 


Cu/C 


450 


0.855 ± 0.035 


0.859 ± 0.009 


W/C 


450 


0.761 ± 0.043 


0.750 ± 0.015 



(b) J/ip suppression ratios pA/pC [22,23] 



Table 1: J/ip suppression in p — ^4 collisions 



We thus conclude that the attenuation of J/ and tp' in p — A coUisions 
is quantitatively well described as the absorption of a pre-resonance charmo- 
nium state in nuclear matter, with aabs = 7.3 ±0.6 mb. The A-independence 
of ip' / {J / ip) establishes that the absorption occurs before the formation of 
physical resonances. 

3. The Centrality Dependence of Nuclear Collisions 

In nucleus- nucleus collisions, charmonium production can be measured as 
function of the centrality of the collision, and hence we have to calculate the 
survival probability at fixed impact parameter h. For an A — B collision, 
it is given by 



d?h AB a(NN -0) 



dajAB V>) 



(fs dz dz' pa{s, z) pB{b - s, z') exp 1-{A-1) J dzA Pa{s, za) c^abs 



X exp |-(-B - 1) y , dzB psib- s, zb) (Jabsj ■ (6) 

Here s specifies the position of the production point in a plane orthogonal to 
the collision axis, while z and z' give the position of this point within nucleus 
A and within nucleus B, respectively. The nuclear density distributions pA 
and pb are defined as above. Note that dSci/d'^b is normalized such that 
when aabs = 0, 

Jd^b j^ dSci{b;aabs = 0) ^^^ 

To obtain a normalized probability at fixed impact parameter b, we have to 
divide [dSGi{b)/d%] by [dSaiib; aabs = 0)/d%]. 

Experimentally, the centrahty of the coUision is determined by a calori- 
metric measurement of the associated transverse energy Et, i.e., the total 
energy emitted in the form of hadrons in a plane orthogonal to the beam 
axis. The more central a collision is, the more nucleons from target and 
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projectile will participate in the interaction, and hence the more secondary 
hadrons will be produced. Starting with the P6-beam experiments at CERN, 
the NA50 collaboration uses in addition a zero degree calorimeter; it counts 
the number of projectile nucleons which have not participated in the collision 
and therefore provides a further check on the centrality of the reaction. To 
compare our calculations with the measured survival probability S'^^{Et), 
we thus have to establish and test a correspondence between the impact pa- 
rameter b, the transverse energy Et, and the energy Ez reaching the zero 
degree calorimeter. 

This correlation is given in terms of the number of "wounded" nucleons 
|25j . Any given nucleon inap — AoiA — B collision is designated as 
wounded when it has interacted inelastically one or more times with one or 
more other nucleons, and fragmentation of such wounded nucleons eventually 
produces the observed secondary hadrons. At high energies, a projectile 
nucleon in general passes the entire target nucleus in a time much shorter 
than needed to transfer information across its own spatial extension, and 
vice versa for the target nucleons. For this reason, the number of times a 
nucleon is wounded does not affect the subsequent production of hadronic 
secondaries in its fragmentation. A given wounded nucleon produces on the 
average Nh hadrons, and if each of these hadrons has an average transverse 
energy qt, then each wounded nucleon will contribute q = Nhqn to the overall 
transverse energy produced in the collision. We thus have the relation 

Erib) = q N^{b) (8) 

between the average number A^^^, of nucleons wounded in a collision at fixed 
impact parameter b and the associated average transverse energy Et pro- 
duced in that collision. In general, q will depend on the incident energy 
in the same way as the average multiplicity. Besides this, in the analysis 
of specific experimental results, q depends on the details of the detector, in 
particular on the rapidity and transverse momentum range in which the pro- 
duced secondaries are measured. Note that relation (jH} is no longer valid 
when Et > ET{b = 0), i.e., when Et is increased by multiplicity fluctuations 
at 6 = 0. 

The average number of wounded nucleons in an A — i? collision at impact 
parameter b is given by 

Ar^«(6) = AJdhTAis}{l-[l-ar,TB{s-b)f} + 
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+ BJdhTB{s-b){l-[l-ar,TA{s)]^}. (9) 



Here cttv — 30 mb denotes the inelastic nucleon-nucleon cross section without 
the diffractive production component, and Ta{s) = J dz pa{z, s) the nuclear 
profile function; the s-integration runs again over a plane orthogonal to the 
collision axis. The two terms in Eq. © thus correspond to the wounding of 
nucleons from nucleus A in its passage of nucleus B and vice versa. 
Integrating Eq. Q over impact parameter yields 



A J d^s {l - [1 - aM^f} + B J d^s a^TAis}]^} , (10) 

since / d'^b T{s — b) = 1 for any s, see Eq. Since the inelastic cross section 
for nucleon-nucleus collisions is given by 

aA = J d's {l - [1 - a^TAis)]^} , (11) 

Eq. ()10|) can be re-written in the familiar form [23] 

j d%N^^{b) = {AaB + BaA). (12) 

Because of fluctuations in the number of wounded nucleons and in the trans- 
verse energy of the secondaries that each wounded nucleon produces, there 
will be corresponding fluctuations in the relation between Ef and b. We as- 
sume the dispersion D in the number of wounded nucleons to be proportional 
to a//V„m 



D\b = « N^^ib), (13) 

with a dimensionless universal physical parameter a found to be about one 
in A — B collisions ,2^. The Et — b correlation function (see Eq. (jH))) can 
then be written as 

. , ' exp|- l^;;;™n ; (14) 

27iq^aN^B{b) I 2q^aN^^{b) J 
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with 



/ 



dEr PAB{ET,b) = 1, 



(15) 



so that it is normahzed at fixed b. 

We now have to check whether the resulting form indeed accounts for the 
observed E'r-distributions and fix the fiuctuation dispersion a. These general 
questions are best answered using the minimum bias data of the NA35/49 
collaboration that will also allow us to later address charmonium pro- 
duction data with as well-determined a form as possible. The minimum bias 
cross section as function of Et is given by 



where Po{b) = [1 — aNTAsib)] is the probability for no interaction. Hence 
[1— Po(&)] cuts off the integral at large b, when interactions become impossible 
(for a hard sphere model, this would be at 6 = Ra + Re)- In Fig. 3, we 
compare Eq. p6p to the measured minimum bias E^ distributions from S — 
Au and Pb — Pb collisions [2^1 • Note that the S — Au data were measured in 
a larger rapidity window than the Pb — Pb data; they were then reanalysed 
to obtain the 5* — Au ii^-r-distribution in the Pb — Pb rapidity window j2Zj ■ 
Setting a = 1 [25J and with g = 1.5 for 5* — Au and g = 1.4 for Pb — 
Pb interactions, Eq. (fT^ is seen to provide an excellent description of the 
measured Ej- distributions. Since the acceptance is the same in both cases, 
the difference in the observed g-values is expected to reflect the difference 
in collision energies; it thus provides us with some estimate of the energy 
dependence. 

The distribution of the transverse energy associated to charmonium pro- 
duction will differ from that just discussed, since the A — B dependence of 
hard processes is determined by the number of collisions, rather than by the 
number of wounded nucleons. When a nucleon, in passing a nucleus, gets 
wounded several times, at present energies this does not change the multiplic- 
ity of secondaries in its subsequent hadronisation; but each of the collisions 
can in principle produce a hard dilepton pair or a heavy quark state. The 
number of collisions in an A — i? interaction at impact parameter b is given 




(16) 



by 



N^'^ib) 



AB TAsib) aN, 



(17) 
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where T^b(6) = / cPs TA{s)TB{b — s) is the nuclear overlap function, normal- 
ized to unity: / (fb TAnib) = 1- Using it, we can define the "conditional" 
i?T-distribution, associated with a hard process such as Drell-Yan or char- 
monium production, 

^ = ^ fd%N^^ib)PMET,b). (18) 
It is seen to be correctly normalized, giving 

<y\B = — I d% N,{b) = AB a^, (19) 
(Jn J 

where is the cross section for the corresponding hard process in an ele- 
mentary A^A^ collision. 

Before we can compare our charmonium survival probabilities at fixed Et 
with data, we thus have to check that also the conditional i?T-distributions 
are correctly reproduced with our Et — b assignment. Since the acceptance 
in the NA38/50 set-up differs from the minimum bias experiments studied 
above, we have to determine once again the mean energy q measured per 
wounded nucleon; the dispersion in A^^^, is retained as determined above. We 
see in Fig. 4 that with g = 0.75 for 5 — ?7 and q = 0.40 for Pb — Pb col- 
lisions the conditional Et distributions are reproduced very well. We have 
included in our distribution the measured J/ip suppression relative to the 
Drell-Yan continuum, since the presently available distributions are domi- 
nated by J/ip rather than Drell-Yan events. At low Et, the data fall below 
the calculated distribution; this is in both cases due to the Et acceptance 
profile of the experiment |28j. The larger difference between the two g- values 
here, as compared to minimum bias, is due to different rapidity coverages of 
the NA38 and NA50 calorimeters. 

As already mentioned, the Pb — Pb experiment provides full collision ge- 
ometry determination by means of a zero degree calorimeter (ZDC), which 
measures at each Et the associated number of projectile spectators - those 
projectile nucleons which reach the ZDC with their full initial energy Em = 
158 GeV. This additional information uniquely identifies the peripherality of 
the collisions leading to the measured charmonium production; qualitatively, 
the NA50 results are quite similar to the minimum bias results of WA80 
j2ni- The measured Et — Ez correlation jSUj in Fig. 5 indicates, e.g., that 
for transverse energy Et — 90 GeV approximately half the projectile nucle- 
ons remain spectators; therefore the impact parameter in this case must be 
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around b ~ Rpb- In general, the number of projectile spectators is A — A^^, 
so that Ez = {A- N^)Ein. Using Et = q = 2q N^, this leads to the 
Et — Ez correlation shown in Fig. 5 (labelled 'Glauber'); it is seen to agree 
very well with the measured correlation. 

4. Nuclear Absorption and Anomalous Suppression 

With the relation between the measured transverse energy Et and the 
impact parameter b of the collision determined, we can now calculate the Et 
dependence of the charmonium survival probability in nuclear matter. It is 
given by 

Sgi{Et) = J d%PAB{ETM 



(20) 

Here the survival probability at fixed impact parameter is given by Eq. © 
and the Et - b distribution by Eq. (HH). Through Eq. Sgi{Et) is 

normalized at fixed Et such that Sgi{Et) = 1 for aabs = 0. With aabs = 
7.3 ± 0.6 mb, from Eq. (^, and with PAB{ET,b) as determined in the last 
section, we have absolute predictions for the nuclear absorption suffered by 
charmonia in nucleus-nucleus collisions. 

We begin with J/ip production and first consider the i^y-integrated sur- 
vival probabilities in O — Cu, O — U and S — U collisions. In this case, the 
experimental survival probabilities are given by = ((T^^/Ai?)/((Tp^/2) 
and thus fully determined in terms of measured cross sections. As seen in 
Table 2, the Glauber predictions fall well within the error range of the data. 

Turning now to the ii^T-dependence, we list in Table 3a the experimental 
survival probabilites fOl 1^ , obtained with reference to the rescaled p — D 
data of NA51. They are shown together with the corresponding Glauber 
results, which again agree well with the data. 

We thus find that the J/V' suppression observed in the collision of O and 5* 
beams with nuclear targets can be accounted for by pre-resonance absorption 
in nuclear matter. For ip' production, this is no longer the case. The data 
for the survival probabilities are shown in Table 3b. They are well below 
the nuclear absorption predictions, and the additional suppression increases 



dSGijb) 



X 



d%PAB{ET,b) 



dSGljb; (Tabs = 0) 

d% 
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nexp 


Sgi 

(7.3±0.6 mb) 


Qcorr 

(7.3±0.6 mb) 


0-Cu 


0.57 ± 0.06 


0.63 ± 0.02 




0-U 


0.53 ± 0.05 


0.52 ± 0.03 




S-U 


0.46 ± 0.05 


0.49 ± 0.03 


0.47 ± 0.03 


Pb-Pb 


0.27 ± 0.02 


0.39 ± 0.03 


0.37 ± 0.03 



Table 2: J/ip suppression 'm. A — B collisions [2,10]. 

with increasing Et- We shall return in Section 5 to the interpretation and 
description of this enhanced ip' suppression; here we only want to note its 
effect on J/ip production. It is known that 8 ± 2 % of the observed J/ip^s 
are due to ip' decay [18 . Since the latter is suppressed in S — U collisions, 
the corresponding fraction of the observed J/ip^s must be suppressed as well. 
We thus obtain as corrected J/ip suppression 

8^/^; = 0.92 Sgi + 0.08 S^r, (21) 

the corresponding values are included in the Tables and are seen to generally 
bring data and nuclear absorption results even closer together. 

As a further check with data of higher precision, we consider the directly 
measured ratio of J/ ip to Drell-Yan production. The ii^T- dependence of this 
ratio is predicted as 

[^f/^DY]iET)=GSr/';{ET), (22) 

with G = a^^/a^Y- Since the nucleon-nucleon ratio at 200 GeV/c beam 
momentum is experimentally not determined with the same precision as the 
nuclear ratios, we fit the form ()22j) to the data piU |M] and check whether 
the resulting Et dependence is correct. The result, with G = 46.0, is shown 
in Table 4a and seen to give excellent agreement with the data. 

Next we turn to the new data for J/ ip production in Pb — Pb collisions. 
The integrated J/ if) survival probability as predicted by Glauber theory, with 
the pre-resonance nuclear absorption cross section is included in Table 
2, together with the corresponding data; data and prediction for the Et- 
dependence are hsted in Tables 4b and 5a. We note that the data decrease 
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< Et> 
GeV 




Sgi{Et) 
(7.3±0.6 mb) 


(7.3±0.6 mb) 


25 


0.54 ± 0.05 


0.58 ± 0.02 


0.56 ± 0.03 


42 


0.49 ± 0.05 


0.52 ± 0.03 


0.49 ± 0.03 


57 


0.45 ± 0.05 


0.48 ± 0.03 


0.46 ± 0.03 


71 


0.43 ± 0.05 


0.45 ± 0.03 


0.43 ± 0.03 


82 


0.42 ± 0.04 


0.43 ± 0.03 


0.41 ± 0.03 



(a) J/il) suppression in S-U collisions 



< Et > 
GeV 


^r(^T) 


Sgi{.Et) 
(7.3±0.6 mb) 


25 


0.35 ± 0.07 


0.58 ± 0.02 


42 


0.25 ± 0.05 


0.52 ± 0.03 


57 


0.19 ± 0.04 


0.48 ± 0.03 


71 


0.13 ± 0.03 


0.45 ± 0.03 


82 


0.09 ± 0.02 


0.43 ± 0.03 



(b) V^' suppression in S-U collisions 



Table 3: Charmonium suppression in S-U collisions [NA38]. 



< Et > 
GeV 




46.0 S^J^iEr) 
(7.3±0.6 mb) 


25 


25.2 ± 0.7 


25.7 ± 1.2 


42 


22.8 ± 0.6 


22.8 ± 1.2 


57 


21.0 ± 0.5 


21.0 ± 1.2 


71 


20.2 ± 0.4 


19.6 ± 1.2 


82 


19.2 ± 0. 1 


18.7 ± 1.2 



{&) S-U collisions [NA38] 



< Et > 
GeV 


[{J/^)/DY],,p{Et) 


46.0 SjY^^iEr) 
(7.3±0.6 mb) 


35 


17.8 ± 2.2 


20.2 ± 1.2 


59 


13.2 ± 1.0 


17.9 ± 1.2 


88 


12.7 ± 0.8 


16.6 ± 1.2 


120 


11.4 ± 0.8 


15.6 ± 1.2 


149 


8.6 ± 0.8 


15.2 ± 1.2 



(b) Pb - Pb collisions [NA50] 



Table 4: (J/'^)/Drell-Yan ratios in nuclear collisions. 
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< Et > 
GeV 




SGiiEx) 
(7.3±0.6 mb) 


(7.3±0.6 mb) 


35 


0.40 ± 0.06 


0.47 ± 0.03 


0.44 ± 0.03 


59 


0.30 ± 0.03 


0.41 ± 0.03 


0.39 ± 0.03 


88 


0.29 ± 0.03 


0.38 ± 0.03 


0.36 ± 0.03 


120 


0.26 ± 0.03 


0.36 ± 0.03 


0.34 ± 0.03 


149 


0.20 ± 0.02 


0.35 ± 0.03 


0.33 ± 0.03 



(a) J/ip suppression in Pb-Pb collisions 



< Et> 
GeV 


fm 


Sgi{Et) 
(7.3±0.6 mb) 


51 


0.12 ± 0.02 


0.43 ± 0.03 


84 


0.11 ± 0.02 


0.38 ± 0.03 


115 


0.08 ± 0.02 


0.36 ± 0.03 


144 


0.06 ± 0.01 


0.35 ± 0.03 



(b) ^p' suppression in Pb-Pb collisions 



Table 5: Charmonium suppression in Pb-Pb collisions [NA50]. 



well below the predictions: while the survival probability measured at the 
lowest Et value is close to pre-resonance absorption, that at the highest 
Et is about 40% smaller. We thus have a clear and -E-r-dependent onset of 
additional J/ip suppression. 

The essential result of this section is thus summarized in Fig. 6: for S — U 
collisions, nuclear absorption fully accounts for the observed attenuation, 
for Pb — Pb collisions there is an additional and in this sense "anomalous" 
suppression. 

For the an enhanced suppression had already been seen in S — U 
collisions. In Pb—Pb interactions, the suppression beyond nuclear absorption 
increases even more, as shown in Table 5b. It also continues to become 
stronger with increasing Et- 

5. Suppression by Hadronic Comovers 

Before we can use charmonium suppression as probe for deconfinement, 
we have to determine to what extent suppression beyond nuclear absorption 
can be accounted for by dissociation in a "normal" confined medium. Con- 
sider a charmonium state formed at the center of the two colliding nuclei at 
the moment of their complete overlap, at rest in the center of mass, i.e., at 
xp = 0. For the SPS c.m.s. energy of 17 - 20 GeV, the nucleon distribution 
within each nucleus is Lorentz-contracted along the beam axis to a disc of 
some 1 - 2 fm thickness. This implies that after about 0.5 - 1.0 fm, most of 
the nucleons have swept over the nascent charmonium state, resulting in the 
nuclear absorption described above. After this time, it finds itself within the 
medium produced by the collision, and it is this medium we want to probe. 

The initial density of hadronic comovers is in general determined by the 
density of wounded nucleons [SH]- Since the integrand [dNu,{b, s)/d'^s] of Eq. 
dni) is the local density of wounded nucleons at fixed impact parameter, the 
average density associated to hard collisions is given by 

iT-wib) = J <fs 

Convolution with PAB{ET,b) will convert expression (j^^ into the corre- 
sponding density at fixed Et', the resulting behaviour is shown in Fig. 7. 

Because of the very low binding energy (~ 60 MeV), the ip' could readily 
be dissociated by hadronic comovers. To study this, we first express the 



{b,s) 



dN^ 



ib,s) 



dh 



dN, 
d^s 



ib,s) 



(23) 
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initial comover density in terms of the wounded nucleon density, 

(24) 

where Nh is the number of comoving hadrons produced by each wounded 
nucleon. We then have to take into account that comover absorption oc- 
curs after the absorption of the pre-resonance state; in other words, only a 
ip' which has survived the pre-resonance passage through the nuclear medium 
can be absorbed by comovers. The survival probability of a ^Z^' in a comover 
medium undergoing isentropic longitudinal expansion therefore becomes 

S'^%b) = I (fs e^Y>{-v^h^'Nh[dN^{h,s)/(fs]\og{Nh[dN^{h,s)/(fs]/nf}} x 

where ahip' denotes the cross section for ip' break-up by collisions with hadronic 
comovers, v the average relative velocity between the ip' and the colliding 
hadrons, and nj the (universal) freeze-out density for comovers. In Eq. ()25p 
we have assumed the cross section for ip' dissociation by hadrons to be energy- 
independent; since 2AfD — M^i ~ 60 MeV, autp' is expected to attain its 
geometric value of about 10 mb very near threshold. Rather than estimate 
the various quantities involved, we prefer to use the ip' suppression m. S — U 
collisions to determine the two constants Ci = Nh/uf and C2 = vah^'Nh em- 
pirically and then verify that the resulting values are of the right magnitude. 
In Fig. 8a we see that with ci = 1.3 fm^ and C2 = 0.8, such a comover de- 
scription indeed agrees quite well with the ip' suppression observed m. S — U 
collisions. The density of wounded nucleons in a nucleon-nucleon collision is 
about 1 fm~^; taking this as nj, and using ~ 1.5, we would get Ci ~ 1.5. 
For a h — ip' cross section of 10 mb, C2 = 0.8 fm~^ implies a hadron velocity 
of about 0.5. All constants thus lead to reasonable values of the involved 
quantities. 

To compare this comover absorption to the (lower energy) Pb — Pb data, 
we rescale Ci and C2 by the factor 1.4/1.5 found above in minimum bias col- 
lisions. In Fig. 8b we see that the result agrees also with the Pb — Pb data. 
Finally, in Fig. 8c, we show the combined S — U and Pb — Pb data as func- 
tion of the average comover density ncoiEx) = Nh nyj{ET) in charmonium 
production; the observed behaviour scales quite well in rico- 



dN^{b,s) 
d^s 
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The extension of such a comover approach to anomalous J / ip suppression 
is confronted with two problems. Heavy quark QCD calculations |32j ex- 
clude J/ break-up by interaction with hadrons in the present energy range, 
since the dissociation cross section is strongly damped near threshold; this 
remains true also when target mass effects are included |33J. However, it is 
not a priori clear that charm quarks are sufficiently heavy to assure the ap- 
plicability of heavy quark theory. Although this theory correctly reproduces 
J/'0-photoproduction data, it can and certainly should be checked directly 
|34j . Until such a direct experimental confirmation is given, it remains of 
interest to assume less or no threshold damping; the J/ip will then interact 
with a hadronic comover medium through much or all of its geometric cross 
section. 

The second problem is that the suppression by hadronic comovers in- 
creases continuously; there cannot be a sudden onset. Hence it is necessary 
to check if the presence of comover absorption in S—U collisions is compatible 
with the constraints provided by all presently available data. The J/ip ab- 
sorption by hadronic comovers is readily obtained from Eq. ()25|) : we just 
replace the ■^/''-hadron break-up cross section by that for J/ ip dissociation. If 
we assume this cross section to be around 1.5 mb (i.e., approximately half its 
geometric value) essentially from threshold on, we obtain as much suppres- 
sion as is observed in Pb — Pb collisions, as seen in Fig. 9a. However, with the 
break-up cross section thus fixed, we now have a parameter-free prediction 
for the S — U suppression, and it is evident from Fig. 9a that this does not 
agree with the data. Fig. 9a also illustrates quite clearly that a precise study 
of normal vs. anomalous behaviour is needed to confirm this conclusion. 

We therefore turn to the ratio of (J/'?/')/(Drell-Yan) production in S — U 
and Pb — Pb collisions (Table 4). In Fig. 9b we see that with the absorption 
cross section fixed at 1.5 mb, the J/ip suppression remains in good agree- 
ment with the Pb — Pb data; however, the corresponding comover absorption 
clearly misses the S — U data. In other words: if we tune the J/ ■i/'-comover 
dissociation cross section such as to obtain the amount of suppression ob- 
served in Pb — Pb collisions, then the only way the data have of showing 
dissent is to disagree with the corresponding S — U prediction; and this is 
what they do. 

Finally we show in Fig. 10 the ratio of J/tp to Drell-Yan production in 
S — U and Pb — Pb collisions as function of ra^^,. One way to relate the Pb — Pb 
data with riw (filled triangles) is obtained by applying our b — Et assignment 
to the measured < Et > points, using the Woods-Saxon nuclear distribution. 
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For a second way (open triangles), we use the B — Et assignment obtained by 
NA50 [HI from the i?T-distribution including the low Et detector deficiency 
and using a hard core nuclear distribution. We note that the precise position 
of the peripheral data points as function of is rather dependent on the 
details of the collision geometry and the nuclear distribution. Hence here 
conclusive results may have to wait for higher statistics and more complete 
low Et Pb - Ph data. 

Let us now comment briefiy on some proposed comover descriptions of 
anomalous J/ ip suppression. The first such study including the Pb — Pb 
data postulates ad hoc a comover density different from that in Glauber 
theory, setting 

ncoiEr) =cEt , (26) 

in terms of the hadronic transverse energy. This relation, although correct in 
the fluctuation regime, is not in accord with the peripherality of the actually 
measured collisions. The comover density as function of Ez is given directly 
by Glauber theory, relating the number of unwounded to that of wounded 
nucleons. Combining this with Eq. (j^Uj) gives the E^ — Ez correlation pro- 
vided by ref. [12]. It is included in Fig. 5 and seen to disagree strongly with 
the measured form. The model thus does not reproduce the actual colli- 
sion conditions of the experiment. It also uses aabs = 4.8 mb, which is not 
compatible with Eq. (jl)). 

A subsequent, more detailed comover investigation [T3] uses Glauber the- 
ory consistently throughout. To account for comover interactions, it requires 
the corresponding cross sections for J / ip-hadron and ■^/''-hadron interactions; 
in addition, it takes into account J/ip —>■ tp' and tp' —>■ J/ip transitions. To- 
gether with the pre-resonance absorption cross section, one thus has a total of 
flve adjustable parameters, to be tuned such as to account for the anomalous 
J/ ip suppression in Pb — Pb collisions as well as for the S — U and p — A data. 
Let us illustrate what this leads to. Our (-i/^'-corrected) nuclear absorption 
results in Table 4a flt the data listed there with a x^/d.f. of about 1. The 
J/ip behaviour predicted in ref. [T3] leads to a /d-f- of about 4 and is thus 
excluded with better than 99% confldence level. In our opinion, ref. ^15j thus 
establishes nicely that even after the introduction of four further parameters, 
a comover picture can reproduce the J / ^z^- suppression observed in Pb — Pb 
collisions only at the expense of giving up an acceptable account of the S — U 
data. 
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We therefore conclude that it is not possible to describe consistently in 
terms of hadronic comovers both the "normal" (pre-resonance) J /if) sup- 
pression observed m S — U collisions and the anomalous J/ip suppression 
in Ph — Pb interactions. In contrast, the measured ip' suppression is well 
described by hadronic comover absorption. 

6. Suppression by Colour Deconfinement 

For the fate of a J/ip in dense matter, the basic difference between co- 
mover absorption and colour deconfinement is that the former is always 
present, changing smoothly with kinematic conditions, whereas the latter 
has no effect below some critical threshold. For colour deconfinement, the 
interaction region thus consists of a "hot" (inner) part, where a J/ip can be 
dissociated, and "cool" (outer) part leaving it intact. In the case of small 
nuclei or very peripheral collisions, no region is hot enough, and so there 
will be only pre-resonance nuclear absorption. Hence the relative fraction of 
hot interaction region becomes the determining variable for deconfinement 
suppression [24] . 

As already noted, the J/ip plays a crucial role in probing colour decon- 
finement by charmonium suppression: while the weakly bound ip' is easily 
broken up in confined as well as in deconfined media, only hard and hence 
deconfined gluons can resolve the structure of the small J/ip and overcome 
its binding energy 2Md — ~ 0.64 GeV [321 ■ We therefore concentrate 
here on J/ip suppression. Since it is not clear from Fig. 9 if the density of 
wounded nucleons (and hence also the initial energy density) is a suitable 
variable to describe the observed J/ijj suppression pattern, we do not want 
to assume an equilibrated system. 

The precursor of a quark-gluon plasma is a system in which the partonic 
constituents are no longer distributed in the way they are in hadrons. A 
change in the gluon distribution can arise by interactions between wounded 
nucleons. The amount of internetting between wounded nucleons is thus cru- 
cial: how many collisions does a wounded nucleon undergo? This is measured 
by the interaction density [HH] 

Nc , X 

^ > 0.5, (27) 

where iVc and A^„, denote the number of nucleon-nucleon collisions and that 
of wounded nucleons, respectively. Since in p — A collisions = Ai'c + 1, k 
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here lies in the range 0.5 < k < 1, with k = 0.5 for p — p interactions. It can 
become larger than unity only in A — B collisions, since there interactions 
between wounded nucleons are possible, and it increases with the centrality 
of the collisions. Since Glauber theory provides both A^^ and N^, k is readily 
calculable. 

The interaction measure k is of particular interest because it can be de- 
termined in a model-independent way from experimental data. The number 
NijY of Drell-Yan pairs produced in an A — 5 interaction is directly propor- 
tional to the total number of nucleon-nucleon collisions, while the number 
of produced hadronic secondaries (or the transverse hadronic energy E^) is 
determined by the number of wounded nucleons. We thus have k ~ Nov/Nh 
in terms of measurable quantities and can cross-correlate charmonium sup- 
pression at a given Et (and hence at a given degree of internetting) with other 
observables which might show a rescattering dependence, such as ratios of 
hadronic secondaries. 

For the onset of deconfinement, we assume that once k reaches some 
critical value Kc, there are enough interactions between wounded nucleons 
to provide the hard gluons necessary for J/ifj dissociation. Consider K{h,s) 
at fixed impact parameter 6, as function of the profile s of the interaction 
region. The condition k(6, Sc) = fixes the extension Sc of the hot volume 
at the given h. We now evaluate the J/ip survival probability (jHl) with the 
constraint that the integrand vanishes for s < Sc- The resulting Sfj'^{b) we 
convolute with P{b, Et) to obtain the predicted survival probability at fixed 
Et- 

Before we can compare this result to the data, we have to take into account 
that approximately 40% of the produced J/ip^s come from Xc decay (see ref. 
[IHj for details), and the Xc, with its larger radius (r^ ~ 0.4 fm > ~ 0.2 
fm) and smaller binding energy {AE~^ ~ 0.3 GeV < AE^ ~ 0.6 GeV), is 
more easily dissociated than the J/ijj. It is not clear if x-hadron interactions 
can be treated by heavy quark QCD in the same way as those for J/ip^s; 
however, the formalism does give reasonable results for the x well 
We thus assume that there are two distinct Kc values: k} and nf > n}. Since 
the S — U data show no anomalous suppression, has to be at least as large 
as the maximum value of K,{b, s) in S — U collisions. From 

= [n{b = 0,s = 0)],^ ^ 2.3 (28) 

we thus get a possible starting point for anomalous suppression. The value 
of is considered as open parameter and determined such as to get best 
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agreement with the data. The overall J/ip survival probability Sj^'^lEr) is 
now obtained by combining 40% ^-suppression with 60% suppression, the 
latter reduced by the suppressed ip' component. In Fig. 11 we see that with 
= 2.9, we do in fact get quite good agreement with the results of both 
Table 4 and 5. 

The natural variable determining the amount of suppression in this pic- 
ture is the fraction fj/^ of J/ip^s inside the hot K-region, compared to the 
overall interaction volume. It is given by 

fj/4b) = l dh N,{b,s) SGi{b,s) / N,{b,s) SGi{b,s) (29) 

The J/ip survival probability then becomes 

sf;4b)/Sai{b) = [i-fj/m; (30) 

again the convolution with P{Et, b) provides the corresponding expressions 
as function of Et- In Fig. 12 we compare the data as function of fj/^ with 
the prediction (jHOI); the agreement is quite good. 

In jTT], a similar approach is presented with the initial energy density 
instead of k as relevant variable. Here also it is the relative amount of hot 
interaction region which determines the suppression pattern. An interesting 
alternative starts deconfinement at the percolation point of transverse string 
areas [T!^ . 

In summary, the observed anomalous J/ ip suppression pattern in Pb — Pb 
collisions agrees quite well with a deconfinement interpretation. 

7. Conclusions 

Here we want to summarize our conclusions and list the further experi- 
ments which we believe are needed to corroborate them. 

• Present p — A data determine the cross section for pre-resonance char- 
monium absorption in nuclear matter as (Jabs — 

7.3 ± 0.6 mb. 

• The J/ip suppression in O — Cu, O — U and S — U collisions is fully 
accounted for by the corresponding nuclear absorption; in contrast, 
Pb — Pb collisions show additional ( "anomalous" ) suppression increas- 
ing with centrality (see Fig. 6). 
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• production is suppressed beyond nuclear absorption already m. S — U 
collisions; this increases further in Ph — Ph collisions. 

• Anomalous J/ ip suppression present only in Pb — Ph collisions rules out 
any explanation in terms of hadronic comovers; all proposed models are 
inconsistent with some of the available data. 

• Deconfinement, with an onset at the interaction density of central S — U 
collisions, is in accord with all data. 

In two steps of the analysis, further data are clearly necessary for definitive 
conclusions. 

• So far, S — U collisions provide normal, Ph — Ph collisions anomalous 
suppression. Data from an intermediate A — A experiment straddling 
the critical divide is needed to check the onset and help in specifying 
the relevant scaling variable. A Ph — Ph experiment at lower energy 
(around 100 GeV/c beam momentum) would provide an excellent cross 
check, if sufficient statistics can be obtained. 

The confirmation of a well-defined onset would establish critical behaviour. 
It is then necessary to confirm that this is colour deconfinement. 

• According to heavy quark QCD calculations, quarkonium dissociation 
is effectively not possible in a confined medium of accessible temper- 
atures. To check this result for charmonium, a direct study of J/ip- 
nucleon interactions is needed ("inverse kinematics" experiment j3lj ) . 

Besides these two really fundamental questions, there are several other points 
where further data could help considerably in underpinning the conclusions 
obtained here. Higher precision p — A data at 450 or 200 GeV/c would 
reduce the error margin in the nuclear absorption cross section and thus 
result in more precise Glauber theory predictions. If these data provide both 
J/ip and ip' cross sections, also the A-independence of the ip'/ {J/i') ratio can 
be confirmed with still greater precision, thus corroborating pre-resonance 
absorption as the underlying mechanism. Finally, a Ph — U experiment could 
extend the range of k (or of the relevant energy density) to higher values and 
thus might provide a first check of a two-stage x ~ '4' suppression pattern. 
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Figure 1: ^^-distributions for nuclear cross section fits: a) relative to pD, 
b) relative to pC . 
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Figure 2: J/%1) survival probabilities in p — ^4 collisions, compared to Glauber 
theory results with (Tabs = 7.3 ± 0.6 mb. 
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Figure 3: Minimum bias Ej- distributions for S — Au and Ph — Ph interac- 
tions pSj, compared to Glauber theory results with q = 1.5 {S — Au) and 
1.4 (Pb-Pb). 
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Figure 4: Hard collision Et distributions foiS — U and Pb — Pb [H] 
interactions, compared to Glauber theory results with g = 0.75 (5 — U) and 
0.40 (Pb-Pb). 



29 



^ 200 
> 

0) 

CD 

I- 
LU 




I I I I \ I I I I \ I I I I \ I I I I 
10000 20000 30000 40000 

EzDc (GeV) 

Figure 5: Et — Ez correlation in Ph — Pb collisions |i30j, compared to the 
Glauber theory correlation. The curve labelled 'G-V is the correlation of 
model [12] ; see section 5. 
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Figure 6: J/tp suppression in S — U and Pb — Pb collisions for (a) the survival 
probability and (b) the ratio of J/ip to Drell-Yan production. The curves 
show the Glauber theory results with aabs — 7.3 ± 0.6 mb. 
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Figure 7: The average density of wounded nucleons in hard interactions for 
S — U and Pb — Pb coUisions. 



32 



1.6 



1.4 



1.2 



0.8 



0.6 



Sexpi^G\+CO 



s-u 



Pb-Pb 



1.6 



1.4 



1.2 



- 0.8 



0.6 



20 50 80 50 80 110 140 

Et (GeV) 



b) 




Figure 8: The experimental survival probability compared to comover 
absorption in S — U (a) and Pb — Pb (b) collisions; (c) shows both data sets 
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Figure 9: The experimental J/'^ survival probability (a) and the ratio of 
J/ to Drell-Yan production (b) compared to comover absorption in Pb — Pb 
and S — U collisions. 
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Figure 10: The ratio of J/ip to Drell-Yan production in S — U (circles) and 
Pb — Pb (triangles) collisions compared to (■^/''-corrected) Glauber theory 
results, as function of ra^; (a) filled triangles use our Glauber impact pa- 
rameter assignment, (b) open triangles that obtained from the experimental 
i?T-distribution |9j. 
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Figure 11: (a) The experimental J/%1) survival probability compared to de- 
confinement suppression in Ph - Pb collisions, with = 2.3 and = 2.9; 
(b) same for the ratio of J/ip to Drell-Yan production. 
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Figure 12: Experimental suppression, compared to (^'-corrected) Glauber 
theory results, as function of the relative fraction fjj^ of hot interaction vol- 
ume. 
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